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Aiming at a description of transport processes where the dynamically generated 
width of the states is potentially large a transport equation beyond the quasiparti- 
cle approximation is derived in first order gradient expansion. An effective particle 
number is identified which is exactly conserved by the coarse grained transport 
equation. Using a test particle ansatz for this conserved quantity allows to rewrite 
the transport equation into equations of motion for test particles. The two-body 
collision terms are formulated in terms of the test particles which gain non-trivial 
renormalization factors due to the coarse graining process. 



1 Motivation 

This italk is based on my recent work §. Due to limited space I refer to this 
workEJ concerning more details and references. 

For the description of heavy-ion collisions the following Boltzmann-type 
equation is commonly used: 

dt + -V, - V,[/(t,f)Vp ) fit,x,p) = /cou (1) 
m J 

where the l.h.s. contains the drift and the Vlasov term and the r.h.s. incor- 
porates the collisions. The central quantity here is the on-shell phase space 
distribution /. "On-shell" simply means that the energy is not an additional 
free variable on which / might depend but is instead given by the appropri- 
ate mass-shell condition. Therefore / depends only on the three-momentum 
(and the space-time variables). I will discuss in a moment why this might be 
an unpleasant restriction for various processes relevant for the description of 
heavy- ion reactions. Before I come to that point I want to discuss briefly how 
eq. (0) can be solved: One usually involves a test particle ansatz 

i 

and deduces from the collisions of these test particles and the evolution of 
their coordinates Xi, pi between collisions. In what follows I will concentrate 
on the evolution part. The equations of motion for the test particle coordinates 
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between collisions are obtained by inserting (Q) in (|l|) and replacing the collision 
term by zero. In this way one obtains the equation 
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This equation can be fulfilled by demanding that both expressions in square 
brackets have to vanish independently for every test particle i. In this way one 
gets the familiar Newtonian equations for the test particle coordinates 

m 

From a practical point of view one expects to get a sufficiently correct expres- 
sion for the phase space density if the number of test particles is reasonably 
large. This traditional transport theoretical framework is capable to describe 
states with energies close to or at the mass-shell point. For heavy-ion col- 
lisions at sufficiently high energy, however, broad resonances come into play 
with decay widths of the order of the typical kinetic energies. In addition, 
also the width due to coUisional broadening might become comparably large. 
In this case the quasi-particle approximation which enters the field theoreti- 
cal derivation of the Boltzmann equation (|^) becomes questionable and one 
has to find extensions of that model. The new central quantity of interest 
will be the Wigner transform (see below) of the two-point function denoted 
by S^{t,x]po,p). Note that this quantity depends explicitly on the energy 
and not only on the three-momentum. It is the purpose of the next section to 
derive a transport equation for without resorting to the quasi-particle ap- 
proximation. Then I will discuss how the obtained transport equation can be 
solved by an appropriate test particle ansatz. One of the questions will be for 
which quantity a test particle ansatz is justified. From the discussion outlined 
above it is already now apparent that in such an extended test particle ansatz 
the energy has to be treated as an independent variable. Thus, this ansatz for 
a so far unspecified quantity has the form 



2 From the KadanofT-Baym equations to the transport equation 

It is obviously appropriate for a talk at this workshop to start out from the 
Kadanoff-Baym equations: 

d Ai 
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D<{1,1') = Jdl [I]'^°'(l,l)i:'<(l,l') + S<(l,l)£'^''(l,l')] (4) 
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0^ + ^) D>{1, 1') ^ Jdl 1) D>{1, 1') + E>(1, 1) Z?-(T, 1')] (5) 



with 

2/) = e(xo - yo) (2;, y) - D< {x, y)] , 
= e(yo-a:o) [i?<(a;,y) -i?>(x,y)] . 

The plus/minus sign refers to bosons and fermions, respectively. Actually 
all information of this non-relativistic system is already contained in the first 
two equations (|^) and (^). In what follows it is more convenient however to 
additionally deal with the equation (||) for the retarded propagator. As usual 
for the derivation of kinetic equations from the underlying field theoretical 
equations I assume now that, concerning the transport phenomena one wants 
to describe, the dependence of all two-point functions and self-energies on their 
respective (macroscopic) center-of-mass variables X — {x+y)/2 is much weaker 
than the dependence on the (microscopic) difference variable u = x — y. In 
this case it is useful to perform a Fourier transformation of all quantities with 
respect to their respective difference variable (Wigner transformation) 

D<{X,p) = jd\e'P''D<{X + u/2,X- u/2) , 

expand in gradients with respect to X and neglect all terms which contain more 
than one derivative in X. I also introduce the generalized Poisson bracket 

[A, B] - dx.AdpoB ~ dp, Adx^B - ^xA VpS + Vpv4 V^B ■ 

By the procedure outlined above the Kadanoff-Baym equations (Q^) turn into 
three (gradient expanded) complex-valued equations. These equations can be 
separated into their real and imaginary parts by introducing the following 
real-valued quantities 

S<{X,p) = ±iD<{X,p) , 
A[X,p) = -2Im5-'(^,P) = ^[D>{X,p) -~ D<{X,p)] . 

From the Kadanoff-Baym equations one gets by gradient expansion the follow- 
ing set of six (real-valued) equations: 

(po-^- RcS''=*) S< = ±it,< ReL»''°* - i[r, S<] + 7[±zE<, A] , (7) 
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PQ--^- ReS"=' \A = T VieD"'"' , (9) 
2m ' 



-t2 

bo-;^-ReE'-°',^] = [r, Re^--^*], (10) 
2 m 

Pa- — - Rct'A RcL>'°' = 1 - iry^, (11) 
2™ / 4 

bo - 1^ - ReS-\ Re^-*] - -i[r, ^] (12) 
2m 4 

where I have introduced the width r = i(S>-i;<) = i(S''<=*-S'*'^). In the end I 
am aiming at decoupled equations for the generalized phase space distribution 
S'^, the spectral function A, and the real part of the retarded propagator 
ReZ)''°*. Obviously one has twice as much equations as quantities which one 
wants to determine. Some of these equations must be redundant. From the 
purely algebraic equations (|^:^) one immediately obtains 

AiX,p) = -, , ffi^) , (13) 

ReD-^iX,p) . Po-&-Ref:-Y^P) ^^^^ 

These expressions automatically solve also eqs. ([lo| , |lp' ). In the quasi-particle 
regime, i.e. for F — > 0, the remaining two equations^,^) turn into the mass- 
shell constraint and the transport equation, respectively. For our case of ar- 
bitrary width F, however, these two equations contain the same information 
provided one realizes by inspection of (|) that the combination ±iS^ — ^S"^ is 
effectively a first order gradient. Since all equations are exact up to (including) 
first order gradients the following replacement on the r.h.s. of (0,||) is allowed: 

[±zE<, ...]«[^5<,...]. (15) 

After this replacement (0) and (^) become identical. Using ( p"3|jl^ ) one finally 
gets the desired transport equation for 5^: 

l-TA [po-^- ReS-\ S<] - I A [F, {po - ^ - ReE'^°')5<] 
2 2m 2 2m 

= TS<Ti^^A{=it>S<Ti^^S>) . (16) 
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Before I will point out how this equation can be solved by a test particle ansatz 
1 seemingly change my topic and discuss for which quantities a test particle 
ansatz makes sense and for which quantities it does not. 



3 Test particles I: A simple example 



Imagine the following situation: One has a theory describing the evolution 
of a quantity f(t,p). This theory, however, is too complicated to solve it 
exactly. (This is presumably easy to imagine.) The only thing one knows 
about that theory is that it conserves a particle number N{t) — Jdp f(t,p). 
By an approximation scheme one obtains from the exact theory the following 
approximate (transport) equation (which resembles a Vlasov equation): 



F{t)^f{t,p)^Q 



(17) 



with a force term F{t) and a renormalization k which both are not further 
specified. Obviously this approximate (transport) equation does not conserve 
the particle number N but instead the quantity N{t) = J dp {1 — K{t, p)) f {t, p) . 
I know discuss two approaches to solve ([l7|), namely test particle ansatze for 
/ and for / = (1 — k)/, respectively: In general, a test particle representation 
for any of the two quantities is given by 



1 1 *^ 

\=-J2S{p-p,{t)) 



(18) 



where L is a normalization constant and M the number of test particles. Ob- 
viously the momentum integral over ( p^ is conserved and given by M/L. 
Therefore, a test particle ansatz for / to solve (|l^) seems to be inappropriate 
since the momentum integral over / (which yields N) is not conserved by (p^. 
Nonetheless, it is instructive to figure out which equations one gets with a test 
particle ansatz for /. In this case one finds 



dn{ t,pi{t)) 
dt 



F{t)-[i-Hc{t,pM)]Mt) 



This yields two equations of motion for one quantity: 



dn dn 



dt dt dpi 



= 



- 



F 
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i.e. an in general overdetermined system. This is the consequence of the fact 
that the integral over /, i.e. N, is not conserved by the transport equation 
(|l7|). Instead a test particle ansatz for / = (1 — k)/ yields 



i.e. one equation of motion 



d 

—5{p-pi{t)) , 



1 - n{t,pi) ' 

Obviously this is a very intricate situation: by construction N is conserved by 
the exact theory but not by the approximate transport equation. Therefore 
a test particle ansatz for / does not make sense. Thus I propose the follow- 
ing strategy: Solve the transport equation by a test particle ansatz for / and 
calculate from the solution the expression dN/dt. Since this expression should 
vanish in the exact theory it provides a test for the accuracy of the approxi- 
mation scheme which has led to the transport equation (|l7|). The lesson from 
this simple example is the following: To solve a transport equation by a test 
particle ansatz one has to find a quantity which is exactly conserved by the 
approximate transport equation. It may appear that this quantity docs not 
coincide with the one which is exactly conserved by the exact (quantum field) 
theory. 

4 Test particles II: Prom the transport equation to the equations 
of motion for test particles 

For the full quantum field theory given by the Kadanoff-Baym equations 
the corresponding particle number is given by 



d^xd 



For appropriately chosen self-energies ( "<i>-derivable" ) this quantity is con- 
served. I will restrict myself to the commonly used two-body collision terms 
(Born terms) which indeed conserve A''. Here e.g. the out-rate is given by 

^""'^^^/(^(^(^^'"^ 5()(p + pi-p.-p3) 

>^livip-P2) ± v{p-P3)f S<{X,p,) S>{X,p2) S>{X,p3) . (19) 
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If the quantity N was also conserved by the transport equation ( |16[ ) a test 
particle ansatz for 5^ would make sense. However, one gets 



-Nit) . - Jd'^Jj^S<K + Jd^.J^- i±^^<S>-^^S<) (20) 
with 

aReS-*(^,p) , Po-^-RcS-'(X,p)c)r(X,p) 



K{X,p) = 



9po r(X,p) 9po 



I have not found any reason why the r.h.s. of (20) should vanish and I strongly 
conjecture that it does not. (Actually it is hard to prove that a quantity is not 
conserved.) This conjecture is supported by the finding that a closely related 
quantity, namely 

d'xd'^p 



N{t)^ I -^S<{t,x-p) 



with S*^ = ^TAS"^ (1 — K) is conserved by the transport equation (16): 
d ^r,^^ f^3^ f d^ f,-^<o> „■v>c<^ 



where the last equation holds for the Born approximation to the self-energies. 
I think it is hard to conceive that both A'' and N are exactly conserved by the 
transport equation. Assuming that N is not conserved a test particle ansatz 
for S'^ is inappropriate. Instead, from the test particle ansatz (^) for which 
is the density corresponding to the conserved effective particle number N one 
deduces the equations of motion for test particles between two collisions: 

= ( dtRet'^' + ^dtT 

1- K \ r 



1-K ^"'^ " r 

-2 _ 

with the off-shellness AEi{t,Xi; Ei,pi) — Ei — ^ - KeY7°^{t,Xi;Ei,pi). The 
obviously new ingredient as compared to the traditional quasi-particle approx- 
imation is this off-shcllness. Its evolution in time can be obtained from the 
equations of motion given above to be 

—AE,^ ^— r. 21 

dt T dt ^ ' 
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One apparent merit of this equation is the fact that test particles automatically 
get back to their mass-shell if the width goes to zero. I note in passing that 
a crucial step to obtain ( |2l| ) consists in the replacement ([l5|). Finally I want 
to comment on the Born terms which in a test particle picture describe the 



collisions between these test particles. Obviously e.g. the out-rate (19) has to 
be rewritten in terms of test particle distributions: 



(27r)4(5(4)(p 




it>iX,p) = 

' d^pi d^p2 d'^Ps 
(27r)4 (27r)4 (27r)4 

with = S'^{X,pj) etc. The fact that S*^ instead of is represented by test 
particles causes rather non-trivial modifications for the collision integrals. In 
turn, disregarding these modifications amounts to solving a transport equation 
which might be rather different from the one which one actually wants to solve. 

5 Test particles III: Comparison of the test particle method and 
Klimontovich's approach to kinetic theory 

Before summarizing I want to comment on the strong formal similarity of the 
test particle ansatz (|^) to the microscopic phase space density 

N 

NM{£,P,t)=J25^^\^-^^it))^^^\p-Ait)) (22) 

i=l 

introduced by Klimontovich i for a classical A^-particle system. I first will 
briefiy outline in which context Klimontovich uses this function: Aiming at 
a statistical description of a classical system of N interacting particles the 
quantity ( p2| ) is introduced as the phase space number density of a specific 
microscopic realization of the macroscopic system. Assuming that the particles 
interact via central two-body forces obtained from the potential ^{\xi — Xj\) 
one obtains Newton's equations for the coordinates of particle i: 

p; = -V,.^$(|f. -x,|). (23) 

Combined with the equation of continuity in phase space dNM / dt = one gets 
dtNM + -V^Nm - V, (dy dV$(|f - x'\)Nm{,S'J, t)VpNM = (24) 
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called the Klimontovich equation. Formally this equation resembles the Vlasov 
equation. However, there is an important difference between the Klimontovich 
and the Vlasov equation: The latter is an approximate equation for the en- 
semble averaged one-particle distribution function / (which can be obtained 
e.g. from (|l|) by neglecting the collision term /cou). In contrast, no statisti- 
cal information has entered ( p^ yet. There, one microscopic realization of 
the iV-particle system is considered. The crucial step which involves statis- 
tics amounts to average (24) over all microscopic realizations which form the 
macroscopic state, i.e. to consider an average weighted by the A'^-particle phase 
space distribution function. In this way one deduces from ( p4[ ) an equation of 
motion for the ensemble averaged quantity Nm which coincides (up to normal- 
ization) with the one-particle phase space density /. Obviously in the resulting 
equation the average over the product of two microscopic phase space densities 
appears (stemming from the two-body interaction part). From ( p^ one can 
also deduce an equation of motion for such a product which in turn involves 
products of three microscopic phase space densities. In this way one generates 
a hierarchy of equations which is analogous to the BBGKY hierarchy of n- 
particle distribution functions. It is the merit of Klimontovich's method that 
the whole hierarchy can be deduced from one equation (|2^) . If in the ensemble 
average of ( p^ ) the average of the product of two microscopic phase space den- 
sities is approximated by the product of averages then one obtains the Vlasov 
equation. In this case all possible correlations are neglected. If such corre- 
lations are retained one can obtain — using a proper approximation scheme 
— e.g. collision terms of Boltzmann type. To summarize, in Klimontovich's 
method the microscopic phase space density (l2^ ) is introduced to obtain (24) 
from which all equations of motion for ensemble averaged products of these 
densities can be deduced. Within proper approximation schemes the well- 
known (ensemble averaged) kinetic equations with and without collision terms 
can be obtained. Now one might turn the procedure around: Starting from 
( p^ ) one might try to find a proper way to solve this equation. This can be 
achieved by the ansatz ( p^ ) which one might be tempted to call "test particle 
ansatz". (This, however, is misleading.) One obtains in this way the equations 
of motion (p3|). Solving these equations e.g. by a numerical simulation yields 
indeed one solution of (|2^). This, however, is not the whole story. One is 
actually interested in the ensemble averaged n-particle distribution functions. 
To obtain these one has to solve (|^ for various initial conditions which belong 
to the same macroscopic state. Calculating averages and fluctuations of the 
various simulations (characterized by the initial conditions) yields the quanti- 
ties of interest. To compare this strategy to the test particle ansatz discussed 
in the previous sections one has to realize that it is the purpose of the latter 
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to yield a practicably applicable method to solve the ensemble averaged Boltz- 
mann equation and not a microscopic equation. Thus one does not start with 
arbitrary initial conditions for the test particle coordinates but instead with 
a large number of test particles such that the initial one-particle phase space 
distribution is reasonably well reproduced by the test particle distribution. A 
second difference concerns the collision terms, i.e. dissipative terms which are 
explicitly present in the Boltzmann equation (|l|) but not in the microscopic 
Klimontovich equation (U) . While the equations of motion for the test parti- 
cles between collisions are — for a classical system — indeed identical to the 
microscopic equations ( p^ ) the test particles in addition scatter on each other 
with the respective scattering cross sections which enter the Boltzmann col- 
lision terms. To summarize, in spite of the strong formal similarity between 
the test particle ansatz (H) and Klimontovich's microscopic phase space den- 
sity ( p^ ) there are important conceptual differences in the two approaches. 
This becomes even more pronounced if genuine quantum corrections like the 
off-shell effects discussed above are included in the test particle approach. 



6 Summary 

1 have presented a scheme to derive from the Kadanoff-Baym equations a 
transport equation beyond the quasi-particle approximation for states with 
potentially large width. 1 have argued that the particle number TV defined 
for the full quantum field theory is not exactly conserved by the transport 
equation. Instead it can be shown that the effective particle number N is 
conserved. I have derived equations of motion for the test particle coordinates 
from a test particle ansatz for S'^. Finally I have presented the modifications 
of the collision integrals once they are used to describe collisions between the 
test particles. Once the transport equation is solved by this ansatz one can 
calculate dN/dt and find out to which amount the conservation of TV is violated. 
This provides a test for the accuracy of the gradient expansion. 
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